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[Translator's note: Neither the references nor the figures referred 

to in this report are included in the translation.] 

1. Introduction /7 * 

The flow conditions around a helicopter blade require the use 
of a model to predict the aerodynamic forces that apply to separated 
states. 

Practical conditions (problems concerning small amplitude motion, 
or large amplitude motion at low frequency) and theoretical ones 
(linearization) have led to the formulation of models in the form 
of sets of differential equations that relate the parameters of 
profile position to the coefficients of aerodynamic stress. 

The first phase of the study consisted of showing the capacity 
of such a model to reproduce the effects of small and large amplitude 
vibrations [1 and 2] . 

In a second phase, it was necessary to find a general model- 
building procedure that, starting with a set of experimental measure- 
ments, covered the profile's entire operating range. This is the 
subject of the present report. 

All along the first two phases, the set of differential equations 
evolved towards a greater ability to describe the phenomena 

*Numbers in the margin indicate pagination in the foreign text. 
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encountered. It is still not completely fixed, and modifications 
that would not seriously undermine the previous work are envisionable . 


The third phase will consist of making a complete aerodynamic 
model of a profile and applying it to cases of real helicopter flight. 

2. Summary of Previous Work [1 and 2] 

The most general form of the model, considering the hypotheses 
made, is the following: 


j 

^ iu < 


. for ^ =1.3 

correspond to coefficients of lift, moment, drag 

= istantaneous incidence at the front quarter rela- ■ 
, ->• 
tive to the instantaneous velocity V 



This model contains a large number of coefficients that have 
been inaccessible to theory for a long time, but whose value can be 
estimated from a series of experiments. Small amplitude sinusoidal 
vibrations are obvious candidates, but other forms of motion can 



also lead to interesting results. 


This model can be simplified for practical applications by 
taking into account certain physical realities and the needs of 
users . 

One particular model was defined on the basis of Kubner's 
coefficients. The system was written in the form: 



Application to the field of helicopter flight has shown the 
necessity of introducing complex poles to take account of certain 
experimental results in the stalled region. The equations were put 
in a slightly different form so as to take account of the different 
conditions that have to be checked. This form is a little distant 

from the general formulation, but the basic principles are still 
respected. 
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This model was applied to the lift and the moment of a profile 
in the OA family. The value of the different coefficients was 
adjusted on tie basis of small amplitude vibration tests in a wind 
tunnel. 

The module obtained was used to reconstruct large amplitude lift 
and moment loops that have compared very well with the experimental 
results. 

In the rest of this report, we will use a slightly different 
form of these equations (see page 7) : 



Apart from changes in the variables, the difference comes from 
the fact that a and r can take on positive values that are inaccess- 
ible to the groupings 2ay and (1 + a^), which permits a broader 
field of application (but even so, a and r must remain positive for 
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The second equation, which takes into account the effects of 
stalling, becomes operative once AC is nonzero, i.e., once the 
incidence reaches the static stalling angle. Experience shovrs that 
stalling does not actually occur until a little later. To take 
account of this fact, one can introduce the function H(6,0), which 
is zero as long as stalling has not really occurred and unity 
thereafter : 


3 


/lO 

Such a correction is perfectly within the spirit of the linear- 
ization that is the basis of the model; If stalling has not actually 
taken place beyond the static stalling angle, it is not possible to 
use the model, which is linearized about a separated flow. 

3. Difficulties in Finding a Model 

The parameters in this type of model are found by using the 
curves of the response to small amplitude vibratory motion as a 
function of frequency. The response of the model is known analytic- 
ally. The parameters have to be chosen in such a way that this 
response is as close as possible to 1he measured one. There is a 
smoothing of the harmonic response as a function of frequency. 

Now, to obtain an equation that is valid throughout a helicopter 
blade's whole operational range, the sets of parameters thus found 
must themselves be smoothed as a function of incidence and then as a 
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function of Mach number. 


This triple complex smoothing (the harmonic response is 
complex: it has an amplitude and a phase relative to iJie original 

motion that created it) is difficult. 

3.1 Choice of Model 


One can find various model forms that reconstruct very well 
the harmonic response of an aerodyneimic coefficient under a 
specific flight condition. We have retained the form in models 
C and D (see preceding Chapter) for the following reasons: 

— The smoothings as a function of frequency are easily 
accomplished because the formulas employed correspond to the 
general behavior of the aerodynamic coefficients. This stems 
from the fact that these formulas come from differential /13 

equations which have been shown to be an approximation of any 
linear system vibrating at small amplitude. The behavior as 
a function of frequency is strongly present, in amplitude and 
in phase. 

— There are not too many parameters, and it is possible 
to separate out the influence of several of them by means of 
certain simplifications that experience allows us. This is a 
consequence of the preceding point: The responses possess 

the experimental behavior, which is covered by a small family 
of parameters. 

— The small number of parameters we use means that each 
of them has a role that is very distinct from the others. When 
a smoothing as a function of frequency at fixed Mach and 
incidence is made, it reconstructs the value of each parameter 
corresponding to the role that it plays. 

Since the natural evolution is known, the parameters themselves 
once again have continuous, and therefore easily smoothable, values. 
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(Even so, pay attention: Stalling for the incidence and the Reynolds 

number for the Mach number can create discontinuities.) 

3.2 Equations and Definition of Variables 

The dot designates the reduced time derivative. We have: 


-h J J tT) ^ ■•‘J B 

ti ^ - 1^ ^ Cq e J M (^,^) 

In real time, the. system becom.es: 

J CjrV / cl ic,'] H(a,e'J 

I Cq, . Cj^ 


with 


A = >*^/[/ 


a' i )//i 

4 '= 

e'- el/// 


(JU = Mach number 

V = wind velocity 

0 ) = real pulsatance 

V = reduced pulsatance, v = ub/v 
t = real time 

T = reduced time, t = vt/b 
b = profile's 1/2 cord 



e = profile's position variable 
e^; appearance of stalling 
0^: mean position during vibration 

0 : amplitude of vibration 

H = function of 6 , 0 taking on the values of 0 or 1 
H = 0 characterizes dynamically attached flow 
H = 1 characterizes dynamically separated flow 
These flow states only correspond to static regions for 0=0. 

= aerodynamic coefficient (see figure 1) 

static value of this coefficient 

static value of this coefficient for an unseparated 
flow, but which has to be defined for the whole 
field covered by the model (even in regions of 
stalling) . 

characterizes stalling /13 

comes from the unstalled aerodynamics 
Cg 2 is the correction for stalling 
C = the mean position during a vibration 
Cg = amplitude of the vibration 

4 , Behavior of the Model 


$): 


J 4- 3 4 - A. sz - 4-C ACj j H (B^ &) 

i C<{. Cf 4 

Conditions respected: 


— At zero frequency. 


one again finds 




AC = C static. 

q q 


— As long as stalling does not occur, the second member of the 
second equation is zero and its solution is C ^2 = 0* Only 
the first equation is active (on the condition that the 
initial conditions for the second equation are indeed 
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Cq 2 (0) = 0 and C ^2 (0) ~ 0/ which is the case if the flow 
remained attached for long enough time that its response had 
the time to dampen . ) 

The function H(0,G) is the only quantity that cannot be chosen /14 
by small amplitude tests. It makes it possible to delay stalling 
until some time after passing through the static stall angle, an 
effect which is observed in all the large experimental loops. 

If an infinitely small amplitude is assumed, tie flow is either 
always separated or always attached, and the equations behave as if 
this function is absent. 

4.1 Role of the Different Parameters in Small Amplitude Responses 

We will illustrate the various roles by a typical example that 
refers to a lift coefficient for a vibration. abotit the front quarter. 

The values of the different coefficients are; 


= 0,103 6 


^ = 0.2 ^ = 0.15 

yi = 0.087 a = 1 

C =-1 


For the 

nonstalled case: 

T = 0.068 

and 

£i^ = 0 

For. the 

stalled case: 

II 

o 

and 

0.180 


?e 

(The reduction of a in the stalled region is generally encountered 
with this type of coefficient.) 

The harmonic responses are traced for the two configurations 
in figures 2 and 6. 


Q 



a. Nonstalled Configuration 


We make the response ^9 * Ca cori^espond to the input 

We have the relationships: 



=- 


{ 


i'- r 

'I 

i 

IB = 


1 

■ 



V 

Examination of the formulas above directly yields the inter- 
pretation of a and s: 


/15 


^ represents the horizontal asymptote for the real part of 
the response ^ /9 (see curve in figure 4) . 

is the slope of the asymptote for the imaginary part of 
(this asymptote passes through the origin) [sic] . See the 
curve in figure 5. 


J modifies the fashion in which the asymptotes are approached: 
small : They are approached rapidly, 

large : They are approached slowly 

(See the curve in figure 3.) 

The formula shows that the 

entire real part is contained between its initial value ^(/d& and 
its asymptotic value a. In addition, the average value ^ 
reached at the reduced frequency 


is 


In the temporal domain, (i characterizes the speed at which the 
steady state is established. The larger 4' is, the more rapid is 
the response to a given level of incidence. 
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Our three parameters ^ are to a large degree uncoupled. 

As in the nonseparated state that interests us here, the experimental 
results are of very high quality. Obtaining them does not pose 
any problem and leads to excellent smoothings as a function of Mach 
number . 

b. Stalled Configuration 


Using the same notation as the preceding section, the response 
is written: 


r 


^ ^ v’.< ^ 

la/ V' 








The parameters are divided into two groups: coming 

from the first equation and (r,a,e) coming from the second. 


/1 6 


These two groups can be treated separately because experience 
shows that the response of the second equation, C -/Q / preponder- 

'Xj ^\j ^ 

ate over the first, Cgj^/6 , in low and medium reduced frequencies. 


'd )l. 


The opposite occurs at high frequencies, where 

-jirr 

Vie deduce from this that for the first group : 


0 while 


— y has a negligible role (see figure 7) 

— <r- and always have their as^miptotical physical sense. 


Notice all the same that the asymptote for the imaginary part 
is reached at a really high frequency (figures 8 and 9) . 
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Role of r and a: 


The first member ^ has a classical response as 

a function of frequency whose minimum amplitude passes through a 
maximum at \/-t- this expression is defined, i.e., if 

a c v^). 

Roughly speaking, r fixes the position of the amplitude peak 
and a its height and width. In fact, however, the second member of 
the second equation also contains v and alters this simple schema. 

But the curves in figures 10 and 11, traced for a realistic applic- 
ation, show that the roles of r and a do take on their classical 
appearance. 

Role of e (figure 12) ; 

The second member of the second equation is the sum of two 
excitations: - -c ^ and -e "d- These expressions were 

obtained after linearization. 

One can see that e makes the overall response increase, but it 
must be mentioned that a reduction in a has the same effect. 

e differs from a by the fact that: 

— a can increase the amplitude peak, but makes it narrower 
at the same time. e does not change its width, but this 
effect is certainly not enough to classify the influence 
of the two variables based on a set of experimental points. 

— The value of e especially determines the phase of the ampli- /II 

tude maximum relative to the profile's motion. The two 

excitations . 4^1 Q and -f))e — 6 are 90° out of 

d& d ^ dB 

phase. The overall excitation therefore has a certain 
phase relative to 0 . This phase keeps reappearing through- 
out the response of the second equation: 
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If e is decreasing: the phase difference tends toward 

-u/2. 

For e = 0: no phase difference. 

For^lt = 'C, I the phase difference is ir/4. 

If e increases: the phase difference tends towards tt/2. 

The models made up to now show that the experimental points 
make a good separation of the effects of a and e possible. 

5. Building a Model 


Four distinct steps are to be considered: 

— modeling the static behavior of 

— modeling in the nonstalled region; finding T 

— modeling in the stalled region; finding 

— verification on large loops. 

5.1 Modeling the Static Behavior 


This concerns finding analytical formulas that reconstruct 
for the whole range of Mach numbers and angles of incidence. In 
particular, one must find an equation for the angle of stall as a 
function of Mach n;imber. 

This phase is more difficult that it appears. To avoid useless 

problems during the model-building process, it is necessary to 

progress with care. The model's harmonic response is going to 

depend strongly on 3C /30/ and one therefore has to approach this 

qs 

term as closely as possible. That does not pose any problems at 

small or very large angles of incidence. Unfortunately, once 

stalling appears, there is a critical region in which varies 

greatly, and therefore 3C /30 equally so. In addition, the 

qs 

experimental data for this region contains large inaccuracies. 

The lack of precision for 3C /30 in the critical region is in /I 

general inevitable, particularly if one is focusing on a moment. 

Even so, it still has to be represented as well as possible. 
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The measurement of C must therefore be conducted very care- 
fully and for a large number of incidences so as to guarantee contin- 
uity in its evolution. A solution consists of using a very slow 
incidence ramp, which yields a continuous measurement. 

The complexity of stalled aerodynamics is such that the use 
of a ramp of this sort is going to give a result that differs 
according to the ramp ' s direction. In the form proposed, tie model 
cannot account for this behavior. The user will have to use a mean 
static curve that seems satisfactory to him. 

5.2 Modeling in the Nonstalled Configuration 

This concerns finding the values of the coefficients t~ of 

Ccj r A = A * (Ayi tt) b -hA & 

throughout the nonstalled region. 

Finding these values is based on the harmonic response of to 
a small amplitude sinusoidal vibration. 

For conventional linearized aerodynamic theories, this response 
is independent of the average profile incidence. Experience approx- 
imately confirms this fact, which induces us to use coefficients 

that depend only on the Mach number. This simplification 
does not always make the equation linear, moreover, since one is 

sometimes led to use a nonaff ine function C 0) C). 

q£ 

We saw in the preceding chapter which roles the different 
parameters play. The quality of the experimental results recorded 
in this region and the very precise physical meaning of each para- 
meter make finding ^ , A, easy to do manually. However, we use the 
algorithm described in Appendix A, which automates this process. 
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5 . 3 Modeling in the Stalled Configuration 


— H(0,0) is a function that we do not know how to choose for 
the moment. A delay depending on the magnitude of 0 when passing 
through the angle of stall 0^ should give good results. In fact, 
we use the universal delay defined by Beddoes [8] for harmonic 
oscillations and v;hich has the value At =; sb/v^ or simply Ax = 5 
in reduced time (exactly 4.88 with a mean deviation of 0.98 in 
142 examples that were chosen by the author and had different pro- 
files, frequencies, vibrational amplitudes, and Mach numbers). 

This delay is the result of the appearance of a vortex on the 
leading edge, which makes the moment drop. The lift does not fall 
immediately, but with a mean delay of 10.8 in reduced time, still 
according to Beddoes, because of this vortex's transit time. However 
this effect is not supposed to be taken into account by the equations 

The choice of the function H(0,0) is independent of the rest 
of the modeling. 

— We also have to determine the values of the six coefficients 
A f ^ , a of the two equations. 

Six parameters is a lot, but by referring to the preceding 
chapter, several simplifications are possible. We had concluded 
that the second equation had a dominant role, except at high frequen- 
cies, where Cg/0 approached the asymptote • -he role of 

in particular was negligible. 

In addition, the first equation must be considered as giving 
C at nonstalled incidences, and even at incidences higher than the 
angle of stall. There exist artifices that delay separation; aspir- 
ation of the boundary layer or high pitching motion or, numerically, 
resolution of the simplified aerodynamic equations. We notice then 
^ a response that is not continuous with that obtained at lower angles 
of incidence. The law = 2ir0 (with 0 in radians) in particular 
seems valid. 



We will therefore try to maintain the first equation's non- 
stalled behavior while only making a few corrections in the para- 
meters . 

We therefore simply propose to: /2Q 

— extrapolate A in the stalled zone. The experimental 

points cannot help to distinguish a particular behavior 
for . As we chose ^(B)^ , the extrapolation is 

all done. 

— slightly correct a and s in the separated region to take 
account of the tendency to high frequency that experiments 
have supposedly shown. 

If the correction does not appear on the experimental points 
(asymptotic behavior reached at too high a frequency, particularly 
for s) , obviously the nonstalled values are retained. 

Finding the Coefficients r, a, e 


In trying to approach as closely as possible the experimental 
points, we define an error function equal to the sum of the squares 
of the errors at each point of measurement. We might at the same 
time assign weights to each frequency, or for the real parts relative 
to the imaginary ones. 

The program in Annex B searches for the values of r, a, e that 
minimize this error for each Mach number and angle of incidence 
(optimization by the deflected gradient method) . 

We could not smooth the r's, a's and e's obtained by functions 
of the Mach number and incidence. However, the smoothings that can 
be done are difficult and lead to poor results because of the dis- 
persion of the optimized values of r, a, and e. 

The action of the three variables must be separated, as described 
in the preceding chapter. 
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The peak ^ /e response is the feature of the most reliable 
curve, measured experimentally. It changes little from one angle 
of incidence to another, and from one Mach number to another. Vie 
have said that its position was mainly defined by r, and equally by 
a, which especially affects its height. 

We will therefore content ourselves at first with smoothing the 
values of r by an expression ^ (Ul,6:) • The new value can be consid- 
ered as a better approximation of the value of r than that obtained 
by gross smoothing. This is the one we V7ill use from now on. 

Still using the program in Appendix B, we then look for the /21 

optimum pair (a,e) , V7ith r fixed at the value that was just deter- 
mined. The different values of a are smoothed by an expression 
& ( &) . This time again, the new value of a can be considered 

as a better approximation of a than those obtained by gross smoothing 
and it will be used from now on. 

All that remains is e, and we also will look for its optimum 
values smoothed as a function of and S 

The model is thus finished. 

5.4 Important Remarks Concerning the Choice of Smoothings 

— We have to find expressions in terms of incidence that 

are valid in the two regions, stalled and nonstalled. These 

expressions must not vary within the nonstalled zone, nor vary by 

very much at high incidences where the flow configurations resemble 

each other. In fact, it seems simpler to use AC (0) as the variable 

in place of 0 . And as certain C variables have too violent a 

y 

behavior, we will rather use AC , which seems the best estimate of 

z 

separation. 

Thus, at AC^ = 0, we find the initial ordinate of the 
curves for ^ J , <3 , e that v/ill be the value of the 

corresponding variable in attached flow. 
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— The functions r and a must never take on negative values. 

That would mean that at a certain fixed incidence, the profile 
would experience forces diverging towards infinity. That would 
be acceptable only in the case of a cylinder, where an oscillating 
force spontaneously develops for certain wind velocities, but it 
would then be necessary to introduce nonlinear terms to prevent this 
force from going to infinity by establishing boundary cycles. (All 
this is developed by Szechenyi [3].) 

— We cannot isolate the behavior of the first equation in the 
separated region, and it was necessary to suppose that it did not 
change significantly. Similarly, the behavior of the second equation 
is not accessible in the nonseparated zone. 

Now, when the flow reattaches, the excitation of the second 
equation cancels out, but not ^^ 2 * second equation continues 

to respond as long as 0^2 is not dampened, which occurs in a manner 
that depends on r and a. 

We therefore also have to define r and a in the attached region, /22 
but small amplitude vibrations cannot inform us in this matter. 

Special experiments are necessary, for example, with levels of 
high amplitude incidence between the stalled and nonstalled regions. 

The manner in which the steady state is established depends on r and a. 

Up to now, we have used large oscillations in incidence center- 
ed on the angle of stall. The fashion in which reattachment occurs 
during the large loops also depends on r and a, but in a more compli- 
cated way. 

A little groping and, of course, the continuity condition for r 
and a at the beginning of separation (but this criterion is not 
very precise) has led us to use up to now: 

r = 0.0100 

a ranging between 0.15 and 0.20. 
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5.5 Summary of the Method 


— Model of static curves 

-- Model for nonstalled incidences (program in Appendix A) 

— Model for stalled incidences 

- finding expressions 

- optimizing triplets (r,a,e), from which is deduced (dL,B) 

- optimizing pairs (a,e), from which is deduced ^ CdL.B) 

- optimizing e, from which is deduced e. 


The optimizations are conducted for each Mach number and angle 
of incidence by the program in Appendix B. 


— Verification in the large loops encroaching on the stalled 
region of the choice of values ^ ( M, and ^ (M, . 

Possible correction of these values with the constraint of remaining 
compatible with the experimental results obtained with small ampli- 
tudes . 

The large loops also make it possible to judge the general qual- 
ity of the model and the necessity of introducing the function H(0,0). 


6. Example of Modeling 


/23 


Example 
motion about 


concerning the lift of the OA 209 profile for pitching 
the front quarter. 
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6.1 Definition of Variables 



8 


9 


P., 


= the profile's incidence in degrees 
= coefficient of lift 

= variables introduced to model the static layers 


The other variables are defined as in chapter 3, page 7. 
6.2 Points of Measurement 


They were obtained by CEAT [Center for Aerodynamic and Thermal 
Studies] [4] . 

The model's field of application goes up to an angle of 
incidence of 20° for Mach numbers ranging between 0. and 0.3, 
but its extension up to Mach 0.4 should be possible. 

Tests with small amplitude vibrations were made at Mach 0.12, 
0.20, and 0.30 for the followiong mean incidences: 0°, 6°, 9°, 

11°, 12°, 13°, 14°, 15°, 16°, 17°, and 18°, and even 20° for 
Mach 0.12. 

The frequencies of vibration were: 

for Mach 0.12: 1.5, 3, 5.5, 10, 15, 20, 30, and 40 Hz 

for Mach 0.20: 3, 5.5, 10, 15, 20, 30, and 40 Hz 
for Mach 0.30: 4, 5.5, 8, 10, 15, 20, 30, and 40 Hz 
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or reduced frequencies ranging between: 0.045 and 1.23 (M = 0.12) 

0.055 and 0.74 (M = 0.20) 
0.050 and 0.49 (M = 0.30) 
Forty Hertz represents a physical limit for the setup that cannot 
be surpassed without harming the quality of the measurements. 


/2 


This set of frequencies made it possible to get a good descrip- 
tion of the harmonic responses necessary for finding the parameters. 


6.3 Static Model 


We looked for the static curves in the form: 






for ^ ^ ^ d' 


-^1 


for d'>^J 


which can be written: 








This curve family is interpreted in this manner (see figure 13) : 

— At low incidences, follows an affine relationship of 

slope Pq. 

— At very high incidences, it tends towards an asymptote of 
slope (and which passes through a point of abscissa 0^ 
and ordinate 

— Finally, in the stalled zone, C tends monotonically towards 
its asymptote at a speed that increases as ^ becomes more 
negative. 
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The values used to reconstruct the experimental curves are: 


Cy = 0.03 
A = 0.102 
T, = 0 

iC =-0.55 M ♦ 0.65 
=-0.3 M -0.43 

= 1.27/ /> ■ = 12.45 \! A- dl'’" 

See the tracings in figure 14. 

A bend in the experimental points is observed at the moment 
when stalling occurs, whereas the smoothings show a sudden stall. 

•We did not take this curvature into account, since that would 
greatly complicate the fromulation to achieve a negligible gain. 

6.4 Model for the ~ Nonstalled Configuration /25 


(The angle of stall is given by the formula dj = 12.45 / A- oU'' ) • 

The program in Appendix 1 enables one to find the set of param- 
eters that makes for the best least squares approximation 

at each Mach number. 

A problem sometimes arises: The optimum curves obtained have 

to pass through the point of ordinate 8C /30 at zero frequency 
in the real part. However, in certain cases inaccuracies cause this 
point to be not properly located on the prolongation of the low 
frequency experimental points. The smoothings are then perturbed. 

In these particular cases, we impose the best point of departure 
possible when applying the program. 

In using the parameters thus determined, we again employ 
3C /30 because the equations absolutely must include the static 
behavior. This shifts the harmonic responses slightly. (See 
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figure 15, for M = 0.12 in particular). 


The retained parameters, which yielded the smoothings in figures 
15, 16, and 17, are the following: 


r/ = 0.20 

^ = 0.087 

^ = 0.0775 -0.08 M 

They will be used for the whole range of incidences lower than 
the angle of stall. Figure 16 shows the effect of incidence on the 
experimental points. This effect, which was neglected, becomes 
important only near the angle of stall, in the zone where the static 
curve has the rounded form that we did not take into account. 

6.5 Model the Stalled Configuration 


The principles that we have established are: 

— A is prolonged by continuity. 

If the experimental points show clear asymptotic behaviors, 
deduce from them the functions artd> 5?^. 

If not, prolong the values obtained for the nonstalled state. 

— Search for r, then a, then e. 

We will therefore use A(ciU,&J= 0.20. 

As the experimental points poorly define the imaginary asymptotic 
direction, we keep = 0.087. 

We now have to choose 'rCc/l,^) . Here, we proceeded in a slightly 
different manner than the established one, v^hich makes the search 
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for the parameter a little more rigorous. However, it is certain 
that a choice based on an estimate of the asymptote of the real 
part would have given a good result. 

a is treated like r, a, and e by the program in Appendix B. 

We have already stressed the fact that the order of finding 
the parameters can. not be random (see page 15) . Since the para- 
meter r has a more important role than a, it must be found first 
in order to stabilize the curves about reliable families. We then 
looked for in order to fix the height of the peaks linked 

to the position of the asymptote and which a and e greatly depend on. 
Then we found a (Ul.ff) and then e 

In fact, we interposed the choice of a between r and a in the 
order r, a, e. 


The results obtained are given in pages 42 to 48 [omitted in 
translator's copy of the original]. 

We expected an acceptable continuity as the Mach number changed, 
but we were forced to note a discontinuity between Mach 0.12 and 
0.20. In contrast, the results at 0.20 and 0.30 are close enough that 
one can use the same formulation for both. 


The difference at Mach 0.12 might arise from a Reynolds number 
effect. Be that as it may, such an effect interferes with the con- 
struction of the model. For the application, we will use one model 
for Mach numbers belov; 0.12, another for those above 0.20, and an 
interpolated model between the two for the intermediate Mach numbers . 


The successive results obtained by the smoothing program are 
given on pages 42, 43, and 44. 


Figure 17 ; We have in fact smoothed /r, which is an approxim- 
ation of the position of the response's peak amplitude when a does 
not take on too high a value. We looked for it among a family of 
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hyperbolas, but parabolas would certainly have been simpler. 


The advantage of hyperbolas comes from their coarser growth 

for very large values of AC , which might be encountered by the 

z 

helicopter blades at certain operating points. In fact, the other 
coefficients were smoothed by parabolas and even a cubic curve, 
which cancelled out the advantage counted upon. For the application, 
one must therefore be careful not to go beyond the smoothed domain. 


with 


= OJ ^ 0( 
/ . - 1 




«' = 1 for 0.12 


K = 1.525 - 4.575 M for 0.12 ^ M ^ 0.20 

0.65 for 0.20 < M^0.40 
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Figure 18 ; The smoothing of a has not posed any problem either. 
A family of straight lines sufficed. These lines' imposed point of 
departure was the nonstalled value of a, which has already been 
determined . 


^ = 0.0775 - 0.08 M + 7 


with 

r = 

- 0,19 

M 

0.12 


r = 

1.3875 ] 

SI - 03565 

0.12<M ^ 0.20 


r = 

-0.079 


0.20CM<0.40 


Examination of the tables clearly shows that the order of the 
choices of r, a, and e cannot be random. In first step, r is easily 
smoothable whereas a takes on erratic values. In the second step, 
meanwhile, it can be smoothed properly. 

Figure 19 ; Very easy smoothing of a by parabolas: 

a = a + 6AC^ 
o z 

The choice of an initial ordinate of around 0.20 is logical, 
but it turns out that a very small variation in this ordinate 
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strongly influences the large loops at the moment when the flow 
reattaches and the response of the second equation is dampened. 

Choosing 0.15, 0.20, or 0.25 does not appreciably change the 
function but can lead to large loops that overlap each 
other. Experience shows that the dampening is in general less /28 

than the critical dampening, which here has the value of 0.20. 

This is the reason for the choice of 0.15 for a (AC = 0). (The 
critical dampening equals a^ = /4r. It is such that any return 
to equilibrium occurs by oscillations for a < a^ and without oscil- 
lations for a > a . ) 

— c 

The curves show the very clear difference in behavior at 
Mach 0.12. 

We choose <3 = 0.15 - 

with S = 1.75 for Mi 0.12 

r = -16.25 M + 3.70 for 0.12-e M $ 0.20 
S= 0.45 for 0.20^ M 0.40 

Figure 20 ; The values of e are very well aligned. 

with <£ = - 2.7 for MC0.12 

= 26.25 M -5.85 f or 0. 1 2 ^ M ^ 0.20 

^ = - 0.6 for 0.20sS M 0.40 

It is remarkable that the value e = 0 is so well suited for 
the initial ordinate, although any value would have been acceptable. 

The model is now completely defined. The smoothings obtained 
in the stalled part are given in figures 21 and 22. The set of 
harmonic responses is shown by figures 23 and 24. Note there the 
discontinuity induced by stalling This discontinuity would have 
been replaced by a sudden and continuous variation if the bend in 
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the static curve had been taken into account. It is remarkable that 
the experimental points clearly show this continuous variation over 
the 2° preceding stalling (see figure 16, page 41 [of the original 
text] ) . 

6 . 6 The Large Loops 

In the following pages, we have traced three series of large 
loops (figures 25, 26, and 27). 

Mach 0.12 frequency: 4 Hz figure 

Mach 0.30 frequency: 4 Hz figure 

Mach 0.30 frequency: 8 Hz figure 

The variation in amplitude is ±6° about a mean incidence of 
between 0° and 20°. The loops about a 20° incidence go beyond the 
framework in which the model was made, but. they do not lead to an 
absurd result. 

The calculations were made both with the rough model and by 
taking the delay in stalling used by Beddoes into account. The 
differences between the two sets of results is rather small, but it 
is certainly in the direction of imporvement. The delay At of 5 
in reduced time is rather small. 

These results should be compared with wind tunnel tests. Such 

a comparison might lead to amodif ication in the parameters a and r 

for AC„ = 0. 
z 

This comparison poses a few problems for us that are outside 
the actual subject of this report. It has therefore been put off 
to Appendix C. 

7. Conclusions 


In this article we have reported on a general method for model- 
ing the response to small amplitude vibrations that is valid for the 



whole domain in which a profile is used. The results obtained from 
an example show a behavior exhibiting a continuity of coefficients 
that validates the choice of equations used. 

Previous work [1 and 2] has shown the possibility of using this 
set of equations for obtaining the behavior at high amplitudes so 
long as we restrict ourselves to relatively slov; movements. Actually, 
these make up the realistic cases. 

We therefore think that our model has entered the application 
phase. 

Work should nevertheless continue by building a model of a pro- /30 
file (0A9) over the whole range of Mach numbers and angles of incid- 
ence encountered by a helicopter blade and for both lift and moment. 
Difficulties in dealing with the moment might appear. They would 
be due to a large dispersion of experimental points. 

Improvements that do not undermine the results obtained up to 
now are still possible. 
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Appendix A; Algorithm for Finding the Coefficients in the 
Nonstalled Zone 


In the nonstalled zone, C = C ; C is the solution of the 

4 4i 4 

euqation: 


/57 


^ 0 f 

■f- A = J 7~ J G 6 


T-F 


then 


e.e, e 




■^z 


ana 






-9 _ 


J 2^ ^ 

F 


d 




■y.y 


Dat's and Meurzec's algorithm [10] is used to optimize the experimen- 
tal points by the ratio: N(v)/D(v). 

Let us note the experimental points 


\ Si] ^ y ^ 


Our family of parameters will better represent the experimental 
points to the extent that one minimizes 


Error = ZT \th _ 


t ■ ' ( X 


where is a positive weighting factor 


ana 


denotes the norm ,/ Xt. ^ 
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We will in fact minimize: 


If one takes a^^ = 1 for a first iteration, and , 

with Dj^ evaluated in the preceding iteration, for the following 
iterations, our minimization series v/ill tend towards the sought- 
after minimiira, if the algorithm is stable. 

We will suppose ^ = Aa f T and continue on. The error becomes 



Error minimization is obtained by making 3 (error) /h A , 
3(error)/3b, and 3(error)/3s all equal zero simultaneously. That 
leads to the following linear system: 





whose solution is reinjected into the denominator 
new iteration. The convergence is almost immediate. 



for a 


/59 


Such an algorithm is generalizable if the variables 
sought for in the form; 


A = Ac (-(Ae) 

L = Ic p 
r = T-c Itc (A &) 


with f./ g., and h. being predefined functions (e.g. 

XX X 

The preceding algorithm is applicable, and the system becomes: 
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2 : 2. 



< a* 


M 


with 




and 



Application of this Algorithm to the Second Equation 


/6G 


The same reasoning leads us to the follovring system: 
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with 


Z n llj = Z 2. V 

* ( 4 / * 


H =* 


1 ■ for the first iteration 

[T-iZm <.'r 


cj and A. 


found in the preceding 


being the values 
.iteration. 


One makes 








then 



This algorithm does not make it possible to take a into 

f\, f\j 

account as we did in the example since the response C /e is no 

Si 

longer in the form of a single ratio but is the sum of two ratios. 
If this sum is reduced to a single fraction, the unknown coeffi- 
cients' products will become involved in a way that is untreatable. 

Gross applicatons of this algorithm have always led to 
divergent calculations here. 
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Appendix B: Optimizing the Parameters by Hinimizing a /61 

Function of Several Variables 


Our problem is to smooth a series of experimental points 

Z-. = X. + jY- by a transfer function (C /eXv)/ which depends on 
K K K 5 

n parameters Yj this case r, a, a, e) . 

To do this, we define an error function; 



is a weighting factor that might be used and that depends 
on the reduced frequency, and y another weighting factor for giving 
more importance to either the real or imaginary parts. 


The method consisting of looking for the roots of the system 


? Srr C ^C) 

7Ti 




n 


leads to systems that are impossible to solve directly. 

We therefore tried to find the minimum error by a gradient 
method built around formulas 22 in page 2.14 of [6]. 


1. Algorithm 


1. Choose a point of departure for i = 1. 

We are acquainted with the systems under study well enough to 
know valid points of departure. 


For a lift, we use: 


= 0.04 

T, = 0.1 


d, = 0.2 

e, = 0 
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2. Analytically calculate the gradient of Err. /62 

3. Calculate the search direction , which is not the one 

with the largest slope (-G^) , but an improved direction that takes 
account of the curvature of the surface z = Err (y , ... , * -his 

curvature is deduced from the previous steps: 



II 


t:., 


If i = 1, one makes T. = -G. . 

1 1 


seven 


Lastly, for a better outcome, 
interactions . 


we reinitialize T. at G. 

1 1 


every 


4. Find i such that the error at the point P. + ET., denoted 
m ^11 

F ) , is at a minimum. 


We therefore look for 


dF(S) 

dZ 


0 , or 


21 (T;) - ( IESl) = 0 


((T^)j represents the jth component of T^.) 

4a. First approximation of 

dEU,)- ^ (Oj lo small 

etc -7 

— (Eq = 10 used) 

Choose the value e 2 such that 1 2/^0 greatest power of 2 less 
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than value will then be divided by 2 so as to have a 

greater degree of safety in the calculation (to avoid going too far 
beyond the sought-after minimum, into a region where Err might 
assume a complex behavior) : 



i»,T( 


INT = integral function 

4 b. Find a limiting interval [^3^ surrounding the solution: /63 

dF 

Starting with the property ^ ( 2 = 0 ) < 0 : 
dP 

— if (E=E ) > 0, divide! by 2 

a! 2 2 

until df/d! < 0; 
dP 

— if ^ (2=!2) > 0/ double 
until dF/dZ > 0 . 

In both cases, save the last two values of Z2 encountered as 
end points to the desired interval [Z^, 

The second type of loop (double Z^...) threatens to never end 

if the minimum Z^ is located in an asymptotic direction. In the 

case of our model, certain optimizations do correspond to an infinite 

value for Z„. There is nothing that can be done about it. 
m 

4 c. Reduce the interval [Z3, Z,^] to [Zg, Eg], by using the 
sign of dF/dZ to choose a smaller interval to surround Zj^. 

4 d. Make a linear interpolation in [Z , Z ] to find Z_, our 

5 6 ' 

best approximation of Z^. 
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5. Calculate the new point "*■ return to 2, 
except if the distance is small enough, in which case the 
covergence is considered to be finished ^ 10~^) . 

2. Results / 6' 

Our test for terminating convergence is very severe 
^ 10 ^) , but there is no risk of divergence near the minimum. 

Therefore, the more difficult the test we employ,- the more we are 
sure of being close to the desired solution. Especially since our 
calculations are made on a machine in the office and are free . . . 


large error 
up to 80 
up to 25 
up to IS 
2 

These results are only indications based on the cases taken up. 


Average number of iterations made: 


small error 


4 degrees of freedom i 20,25 
3 ” '' " i10.12 

2 ■ ” .... ^ 8.10 
1 " ” " 2 


3. Possible Improvements 

The formulas used in 3 significantly improve the outcome of the 
system. Other formulas, proposing better search directions, can be 
used (see in particular [7] ) . 


Note that on page 2.15 in [6] are also found formulas 23, which 
have been used too, but which have given poor results. (There is 
probably an error involved.) 
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Appendix C; Comparison of the Model with Real Life for the Example /6E 
Used as an Illustration 


The tracings of the calculated large loops (figures 25, 26, 27) 
reveal a disagreement with the measured large loops (figures 30, 

31, 32, dotted lines) . A similar result was obtained by Miss 
Falchero [9] . 

1. The Problem of Quality in Wind Tunnel Tests 

The main cause of disagreement seems to come from the fact that 
the large amplitude tests and the small amplitude tests, which come 
from two different series, were made vrith two different mockups of 
che 0A9 profile. 

The first series [5] concerned the measurement of the profile's 
response to various harmonic motions of large, but also small, 
amplitude (at 20 Hz, 30 Hz, 40 Hz). The mockup was then destroyed 
in order to. recover the pressure sensors. Since the modeling process 
requires many more points of frequency, a second series [4] took 
place, which led to the model in chapter 6. 

Naturally, a certain dispersion appeared between the two series 
each time the same quantities were measured: 

— Figure 28: The static curves are a little different, with 

stalling appearing a little earlier in the second series and a lower 

stalled C . A third static curve, obtained at Modane, is also 
z 

represented. It gives still another result, with a particularly 
high angle of stall. 

— The small amplitude vibration measurements common to both 
series (20, 30 and 40 Hz) also show variations. One can see in 
particular, thanks to the tables in figure 20, that the maximum of 
the real part, which was located at between 15 and 20 Hz during the 
second series, was around 30 Hz in the first series, which corres- 
ponds to slightly higher values of r. 
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In a general manner, the wind tunnel tests tend to be dispersed /66 
at the moment of stalling and in the stalled region. Any little 
thing is enough to change, for example, the angle of stall. This 
little thing can be an irregularity in the profile or a three- 
dimensional effect. It is possible that small vibrations are more 
reliable measurements than static curves, but this is not obvious. 


All this is rather annoying since the differences noted are 
significant in the model. The use of the first series' static 
curves coupled with the coefficients deduced from the second series 
greatly improves our loops. (See figures 30, 31, and 32.) 



J- 




The model as we have defined it correctly schematizes the 
behavior of a certain mockup in its wind tunnel. The defects in the 
profile of the wind tunnel and of the processing of measurements 
(customarily called the wind tunnel's "age") are blindly taken into 
account in our equations, which account for the profile's large 
amplitude behavior in the same measurement conditions. These con- 
ditions are chosen to be as close as possible to reality. However, 
wind tunnels are hard to master. Obviously, these problems are those 
of any model based on wind tunnel tests. 


2. Application of the Model to the Large Loops in the First CEAT 
(Center for Thermal and Aerodynamic Studies] Series 


For this application, it is necessary to take account of certain 
results obtained in the first series, at least of the static curves 
that the model is supposed in particular to predict. 

We therefore tried to find formulas for the static curves 
that are similar to those already used: 
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V 


I 




with 


^= 0.05 

T, ^ 0,102 / //. (M>^ 
T,^ 0 


A= 4/3 M + 0.7 


- 0.5 

= 16° - 10 M 


8&£>j 


It is this single correction that yielded figures 30, 31, and 
32. The results are nov7 acceptable and certainly would have been 
even a little better if all the tests had been done in succession, 
especially at very high angles of incidence. 


We observed the effect of the corrections more completely on 
a loop straddling the angle of stall. Refer to figure 33. 


The conditions are: Incidence: 11° ±6 

Mach number: 0.30 

Frequency: 4 Hz 

First loop : Rough application of the model obtained from the 

second series. The result was not good. 


Second loop : Introduction of Beddoes' delay, which prolonged 

the right part of the loop, but did not make an overall improvement. 

Third loop : Use of static from the first series. The loop 

was now a good one. 


Fourth loop : Modification of r to take account of the obser- 

vation gleaned from the tables in figure 29: 


40 



OMERA 


DIRECTION DBS STRUCTURES 

Rapport Technique n° 12/1841 RY 080R 
Juillet 1981 


Algorithme de lissages successifs pour la recherche 
du tnodele semi-empirique de prevision des forces ae- 
rodynamiques ins tationnaires developpe a I'ONERA 





OFFICE NATIONAL D’ETUDES ET DE RECHERCHES AEROSPATIALES 

29, avenue de la Division Leclerc, 92320 CHATILLON DP/MLD 

Direction Scientifique de la Resistance des Structures 


170O/0R 



n 


n 


Rapport Technique n° 12/1841 RY 080R 
Juillet 1981 


Algorithme de lissages successifs pour la recherche 
du mod&le semi-empirique de prevision des forces ae- 
rodynamiques instationnaires developpe a I'ONERA 


L 

L 


J 


J 


H. Loiseau 


Auteur : 


Chef de la Division D. Petot 

"Essais Speciaux" 


Vu et approuve, 


Le Directeur Scientifique 
de la Resistance des Structures 



G. Coupry 



.200091 


FICHE D'IDENTIFICATION 


ORGANISME EMETTEUR ; 



OFFICE NATIONAL D'ETUDES 
ET DE 

RECHERCHES AEROSPATIALES 

29, Avenue de la Division Leclerc, 
92320 Chatillon 


Direction : des Structures 


CLASSIFICATION 


Secret militaire 


Titre : NC 
Fiche : NC 
Document : NC 


Secret industriel 


Non Confidentiel Industrie 


Numero de la fiche programme 


R. 212 A 


Service de I'Etat 
charg4 de I'exdcution du contrat 


STPA 


Num6ro du contrat 


Rapport de fin d'dtude 


Convention n° 80.95.005 
Lot 3 - Poste 3 . 1 


TITRE 


Algorithme de lissages successifs pour la recherche du mod&le semi-empirique 
de prevision des forces aerodynamiques instationnaires developpd a I'ONERA 


AUTEUR (S) : 

(personne physique) 


D. Petot 


DATE : 

NUMERO D’ORIGINE du document 

NOMBRE 


Pages 

Figures 


Juillet 1981 

R.T. n° 12/1841 RY 080R 

73 


10 


RESUME D'AUTEUR 


Ce rapport est relatif aux travaux effectues dans le cadre de la fiche 
R 212A du programme 1980 : mdthodes d'aerodynamique instationnaire adoptees aux 
dtudes d'aeroelasticite des helicopt&res . II concerne la mise au point et le de- 
veloppement d'un mod&le semi-empirique des forces adrodynamiques instationnaires 
agissant sur profil oscillant h. grande incidence. L' algorithme de lissages suc- 
cessifs developpd dans ce rapport conduit d'une fa^on tr&s satis faisante aux 
coefficients du mod&le. 

Ce fait confirme sa valeur et le rend maintenant operationnel . 


NOTIONS D'INDEXAGE 


- Aeroelasticite des helicopt&res - 

- Forces aerodynamiques instationnaires en regime ddcrochd 


Voir liste de diffusion au verso. 


74-20-SA 














O.N.E.R.A. 


RT 12/1841 RY 080R 


f^ge 5 


TABLE DES MATIERES 


n° pages 

. Titre et auteurs 1 

. Fiche d ' identification 3 

1 . INTRODUCTION 7 


2. RESUME DES TRAVAUX PRECEDENTS 7 


3. DIFFICULTES DANS LA RECHERCHE DE MODELE 10 

3.1 Choix du modele 10 

3.2 Equations et definition des variables 11 

4. COMPORTEMENT DU MODELE 12 bis 

4.1 Role des differents parametres sur les rSponses 

a faible amplitude 13 

5. CONDUITE DE LA MODELISATION 16 

5.1 Modilisation du comportement statique 16 

5.2 Modelisation en configuration non ddcrochie .... 17 

5.3 Moddlisation en configuration decrochde 18 

5.4 Remarques importantes concercant le choix des 

fonctions 20 

5.5 Resume de la methode 21 

6. EXEMPLE DE MODELISATION 22 

6.1 Definition des variables 22 

6.2 Points de mesure 22 

6.3 Modelisation statique 23 

6.4 Moddlisation en configuration non dgcrochie .... 24 

6.5 Ifodglisation en configuration dgcrochge 24 

6.6 Les grandes boucles 28 






O.N.E.R.A. 


RT 12/1841 RY 080R 


Page 6 


7. CONCLUSIONS 28 

REFERENCES 30 

A : Algorithme de recherche des coefficients dans le 

domaine non ddcrochd 56 


ANNEXE B : Optimisation des pararaitres par la minimisation 
d’une fonction de plusieurs variables 


ANNEXE C : Comparaison du modele avec 1* experience pour 

I'exemple qui a servi d' illustration 64 








O.N.E.R.A. 


RT 12/1841 RY 080R 


Page 7 


1 . INTRODUCTION 

Les conditions d' ecoulement sur une pale d'helicoptere necessitent 
1 'utilisation d'un modele de provision des forces aerodynamiques s'appliquant 
aux regimes decolles. 

Des conditions pratiques (les problemes concement des mouvements h. 
petite amplitude, ou des mouvements ^ grande amplitude, mais de faible fre- 
quence), et thdoriques (linearisation) ont amene U la formulation de modSles 
sous foime de jeux d' equations differentielles reliant les parametres de 
position du profil aux coefficients d' efforts aerodynamiques. 

La premiere etape de 1 'etude a consiste a montrer la capacite d'un tel 
modele a reproduire les effets des vibrations a petite et S grande amplitude 
[1 et 2]. 

Dans un deuxieme temps, il fallait trouver un precede general pour 
I'obtention du modele sur toute la gamme d' utilisation du profil, et ce, & 
partir d'un jeu de mesures experimentales . C'est I'objet de ce rapport. 

Tout au long de ces deux premieres etapes, le jeu d' equation differen- 
tielle a dvolue vers une plus grande aptitude a ddcrire les phenomenes rencon- 
tres. II n'est pas complbtement fige, et des modifications s ont encore envi- 
sageables qui ne remettraient pas en cause les travaux precedents. 

La troisieme etape consistera 5 modeiiser completement un profil et a 
appliquer cette aerodynamique & des cas de vol rdel d'helicoptere. 


2. RESUME DES TRAVAUX PRECEDENTS [1 et 2] 


La forme la plus generale sous laquelle peut se mettre le modele, compte 
tenu des hypotheses retenues, est la suivante : 



3 

+ H 4 * - - - 



= 

pour 

7 = 1,3 


(i> 

^ t 9 

<Pt , 

correspondent aux coefficients de portance, moment, trainee 


t/ = 

incidence instantanee au quart avant par rapport h la vitesse 
instantanee M 
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& s: incidence par rapport a la vitesse moyenne V, du profil 


y* ~ parametre de vitesse sans dimension 

/iT/ 


Les sont des coefficients, 

O 


Modele A 


Ce module comporte un grand nombre de coefficients inaoKssibles S la 
thdorie pour longtemps, mais dont on peut estimer la valeur & partir d’une 
sdrie d' experiences. Les vibrations sinusoidales a petite amplitude s'imposent 
a priori, mais d'autres formes de mouvement peuvent aussi conduire & des 
resultats intdressants. 


Pour les applications, on peut simplifier ce modele en tenant compte 
de certaines realitds physiques et des besoins des utilisateurs . 


Un modele en particulier a ete defini a partir des coefficients de 
K^^ner. Le systdme a dte mis sous la forme : 


nr 


1 + ■) 

if 

\ . 

+ 


\nj 


0 LI 







0 

0 L 






e 


7 




e 


F el M ddsignent la portance et le mouvement 

^ est le deplacement sans dimension le long de la normale 

& est la rotation du profil par rapport au quart avant 

Modele B 


L' application au domaine de vol de I'helicoptere a montrd la ndcessite 
d'introduire des poles complexes pour tenir compte de certains resultats 
experimentaux en zone decrochee. Les equations ont dte mises sous une forme 
legerement diffdrente afin de tenir compte des differentes conditions a 
verifier. La forme s' est un peu dloignee de la formulation generale mais les 
principes de base sont toujours respectes. 




Ce modele a 6td appliqud a la portance et au moment d'un profil de la 
famine QA. La valeur des differents coefficients a §te ajustSe a partir 
d'essais de vibration a petite amplitude en soufflerie. 

Le modele trouve a ete utilise pour retrouver des boucles de portance et 
de moment de grande amplitude et a donnd lieu a de tres bonnes comparaisons 
avec 1' experience. 

Dans la suite de ce rapport, nous utiliseront une forme un peu diffd- 
rente de ces Equations (voir page 11 ) : 



Mis S part les changements de variables, la difference tient en ce que 
i. et A peuvent prendre des valeurs positives inacessibles aux groupements 
z^r et ) , ce qui permet un champ d' application plus libre ; 

(mais a et -t- doivent quand meme rester positifs pour des questions de 
stabilite de la rdponse) . 

La deuxieme equation qui prend en compte les effets de ddcrochage 
intervient des que est non nul, c*est-a-dire des 1' incidence de ddcro- 
chage statique. L’ experience montre que le decrochage n’ intervient effective - 
ment qu*un peu plus tard. Pour tenir compte de ce fait, on peut introduire 
une fonction H ( &!$) , nulle tant que le ddcrochage ne s'est pas vraim.ent 
effectue, unite ensuite : 
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Une telle correction est parfaitement dans 1' esprit de la linearisa- 
tion qui est a la base du modele : au del§ de 1’ angle de ddcrochage statique, 
si le decrochage n'a pas vraiment eu lieu, il ne £aut pas utiliser le modele 
linearise autour d'un ecoulement ddcolld. 


5. DIFFICULTES DANS LA RECHERCHE DE M3DELE 

La recherche des parametres dans ce genre de modele se fait au moyen 
des courbes de rdponse ^ un mouvement vibratoire de petite an^jlitude en 
function de la frequence. La reponse du modele est connue analy tiquement . 

Les paramitres doivent etre choisis de fagon que cette reponse soit aussi 
proche que possible de la reponse mesuree. II y a lissage de la reponse harmo- 
nique en function de la frequence. 

Or, pour obtenir une Equation valable dans tout le domaine de fonction- 
nement de la pale d'helicoptere, les jeux de parametres ainsi trouves doivent 
eux-memes se lisser en function de 1* incidence, puis en function du Mach. 

Ce triple lissage complexe (la reponse harmonique est complexe : elle 
a une amplitude et une phase par rapport au mouvement original qui I'a creee) 
est difficile. 

5.1 Choix du module 

On peut trouver diffdrentes formes de meddles qui restituent trds bien 
la reponse harmonique d'un coefficient adrodynamique dans une certaine condi- 
tion de vol. Nous avons retenu la forme des meddles C et D (voir chapitre 
precddent) pour les raisons suivantes : 

. les lissages en function de la frequence se font aisdment, parce que 
les formules employees correspondent au comportement general des coef- 
ficients adrodynamiques . Cela tient au fait que ces formules proviennent 
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d' Equations differentielles dont on a montre qu'elles etaient une 
approximation de tout systeme vibrant, lineaire a faible amplitude. Les 
comportements en frequence sont bien presents, en amplitude et en 
phase, 

. il n'y a pas trop de parametres, et moyennant certaines simplifications 
que 1* experience nous permet, il est possible de s^parer 1' influence 
de plusieurs d'entre eux. C'est une consequence du point precedent : 
les rdponses possedant le conportement experimental, celui-ci est 
convert par une famille de peu de parametres, 

. le peu de parametres que nous utilisons fait que chacun de ces demiers 
a un role bien distinct de celui des autres. Quand un lissage en func- 
tion de la frequence a ^fediet incidence fixes est effectue, celui-ci 
restitue la valeur de chaque parametre correspondant au role qu'il 
joue. 

L' evolution naturelle etant continue, les parametres se retrouvent eux- 
mones avec des valeurs continues, done aisement lissables ; (attention quand 
meme : le ddcrochage pour 1’ incidence, et 1' influence du nombre de Reynolds 
pour le Mach peuvent erder des discontinuites) . 


3.2 Equations et definition des variables 


Le point ddsigne la derivation par rapport au tenps reduit . On a 

J ^ i-T) ^ ■•■ji S 
Cq^4 ^ A Cef^ = - ^ A. /i Cq -t e ACf J H 
Cf. Cq^*Cq^ 

En temps rdel, le systeme devient : 

J y Cq^ 4-/L ^ *■& A CfJ H (B,6'J 

I 4 ; ■> 


J//J 


a'-- J M 

e'. 
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= nanbre de Mach 
= vitesse du vent 

“ pulsation reelle 
= pulsation reduite ; 

= temps r§el 

= temps reduit ; /l 

^corde du profil 
= variable de position du profil 

Oj : apparition du decrochage 
^ : position moyenne pendant une vibration 
6 : amplitude de la vibration 


H = fonction de $ yd prenant les valeurs 0 ou 1 


H = 0 caracterise I'ecoulement dynamiquement accroche 
H = d •' ” " decroche 

Ces etats de I'ecoulement ne correspondent aux 
zones statiques que pour & =0. 




= coefficient aerodynamique (voir figure 1) 


CiJM. e) 


: valeur statique de ce coefficient 



valeur statique de ce coefficient pour un 
dcoulement non decolle, mais que I’on doit 
definir sur tout le domaine modelise 


(meme dans les zones decrochees) . 
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caracterise le decrochage 
^ provient de 1 ' aerodynamique non decrochee 
^ est la correction due au dScrochage 
position moyenne pendant une vibration 
Cj . amplitude de la vibration 


4. COMPORTEMENT DU MODELE 



J (I * r) & d 

■#- d ^A, ^-c J f! 


Conditions respectces : 

. a frequence nulle, on retrouve 


AC^ = Cf 


, tant que I'on ne decroche pas, le deuxieme membre de la deuxieme equation 
est nul et sa solution est = 0 ; seule la premiere Equation joue (& 

condition que les conditions initiales pour la deuxiane Equation soient 

- # 

bien et Cf^Co)^0 , ce qui est le cas si I'dcoulement est 

reste accrochd suffisamment longtemps pour que sa reponse ait eu le temps 
de s ' amortir) . 
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La fonction est la seule quantite qui ne peut pas se choisir 

par des essais de faible amplitude. Elle permet de retarder le decrochage 
quelque temps apres le passage de 1' angle statique de decrochage, effet qui 
s 'observe sur toutes les grandes boucles expdrimentales . 

Si on suppose une amplitude infiniment faible, I'ecoulement est soit 
toujours decolld, soit toujours accrochd, et les Equations se comportent 
comme si cette fonction dtait absente. 


4.1 Role des differents paramStres sur les reponses faible amplitude 

II sera illustre par un exemple type qui se rapporte a un coefficient 
de portance pour une vibration autour du quart avant. 

Les valeurs des diffdrents coefficients sont : 

Cj( = 0,103 6 /= 0.2 ^ = 0.15 

= 0.087 a = 1 
e =-l 


Pour le cas non decroche : •7“ = 0.068 et = 0 

Pour le cas decroche : 9" = 0 et 0.180 

(la diminution de dans la zone decrochee est en general rencontrde pour 
ce type de coefficient) . 

Les reponses harmoniques sont tracdes dans les deux configurations 
aux figures 2 et 6. 


a. configuration non ddcrochde 
A I'entree 

c,. q,- Q 

On a les relations : 


on fait correspondre la reponse 



iHU-rr- 

Jvr 
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L'examen des formules ci-dessus dome directement 1 ' interpretation 
de r et >4 : 

<T reprdsente 1' asymptote horizontale de la partie reelle de la 

•s# 

rSponse ^ (voir courbe figure 4) , 

S' 

est la pente de 1' asymptote de la partie imaginaire de 
(cette asyn 5 )tote passe par I'origine). Voir courbe figure 5. 

J modifie la fagon dont sont atteintes les asymptotes 

. faible : elles sont atteintes rapidement 

. fl fort : elles sont atteintes lentement 

(voir courbe figure 3) 

La formule montre que la 

partie reelle est toute enti§re comprise entre sa valeur initiale 
et sa valeur asymptotique v . De plus, la valeur moyeme ^ 
est atteinte § la frequence rdduite ))= ri . 

Dans le domaine temporel, J caractgrise la vitesse & laquelle le 
rdgime permanent s'dtablit. La rdponse a un Echelon d' incidence est d’autant 
plus rapide que / est grand. 

Nos trois parametres A , ^ ^ sont ainsi largement decouplds. Corame 

de plus dans le regime non ddcolld qui nous intdrese ici, les resultats exp6- 
rimentaux sont de tres borne qualitd, leur recherche ne pose aucun probleme 
et conduit a des lissages en fonction du Mach, excellents. 



b. configuration dScrochde 


Avec les memes notations qu'au paragraphe precedent, la reponse s’ecrit: 


(S'/ 




V 
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Les param§tres se part agent en deux groupes : issus de la 

premiere Equation et C -t , 2 , ^ ) issus de la seconde. 

On peut trailer ces deux groupes separement parce que 1' experience 

montre que la rdponse de la deuxieme equation est preponderante devant 

celle de la premiere , dans les basses et moyennes frequences 

reduites. Le contraire se produit ^ frequence eievee oil ^ 

tandis que ^/s' 7 ~-t- 

f o jf_, ^ y 

Onen deduitpour le_premier_groupe _( s!" • 

. f/ a un role ndgligeable (voir figure 7) 

. 0 " et ont tou jours leur sens physique asymptotique . 

On remarque tout de meme que 1' asymptote sur la partie imaginaire 
est atteinte a frequence vraiment eievee -figures 8 et 9-. 

Role de _>L_ et_ de_ f _ __ 

Le premier membre + ^ possede une rdponse clas- 

sique en fonction de la frequence, dont 1' amplitude passe par un maximum 
pour (lorsque toutefois cette expression est ddfinie, 

c'est-a-dire si 2 CV^). 


^/b 


B )), 




''En gros", fixe la position du pic d' amplitude et a sa hauteur 
et sa largeur. En fait, le deuxieme membre de la deuxieme equation contient 
)) egalement et vient modifier ce schema simple. Mais les courbes, figures 
10 et 11, tracees dans un cas d'application rdaliste, montrent que les roles 
de -t et de a ont bien leur allure classique. 


Role de _ e _ ^figure 12) 


tions : - <. 

linearise. 


Le deuxieme membre de la deuxieme equation est la somme de deux excita- 

e et -e ^ ^ 
dB 


expressions obtenues apr§s avoir 


On voit que e fait globalement croitre la rdponse ; mais il faut 
dire qu'une diminution de a a le meme effet. 

e differe de a par le fait que 

. a peut faire augmenter le pic d’ amplitude, mais en le 
rendant plus etroit. e ne modifie pas sa largeur, mais cet 
effet n'est certainement pas suffisant pour classer 1* influ- 
ence des deux variables a partir d'un jeu de points experi- 


ment aux. 
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. surtout, la valeur de e determine la phase du maximum 
d' amplitude par rapport au mouvement du profil. En effet, les 
deux excitations --c 2^5 Q et ^ 

en quadrature. L' excitation globale admet done une certaine 
phase par rapport a 9 , phase qui se retrouve dans toute la 
repons e de la deuxieme equation : 

si e ddcroit : le ddphasage tend vers 

pour e = 0 ; pas de dephasage 

pour jJe = ^ : dephasage de JT 

4- 

si e croit : le dephasage tend vers H- 

Les modelisations effectuees jusqu'ici montrent que les points experi- 
mentaux permettent de bien separer les e££ets de a et de e . 


5. CONDUITE DE LA MODELISATION 

Quatre etapes bien distinctes sont a considdrer : 

. model isat ion du comportement statique de 
. moddlisation en zone non decrochde : recherche de ^ 

. modelisation en zone ddcrochee : recherche de <t,A, t, 

. veri£ication sur grandes boucles. 

5.1 Moddlisation du comportement statique 

II s'agit de trouver des £orraules analytiques qui restituent sur 

toute la gamme de Mach et d’ incidence. II £aut en particulier trouver une loi 
donnant 1* angle de ddcrochage en £onction du Mach. 

Cette phase est plus di££icile qu’elle ne le parait. Pour dviter des 
problemes inutiles dans la poursuite de la moddlisation, il est necessaire 
de la conduire avec soin. La rdponse harmonique du modele va ddpendre £orte- 
ment de /^e » et il va done £alloir approcher ce terme au mieux. 

Cela ne pose pas de probleme aux £aibles incidences, ni aux tres £ortes inci- 
dences. Il y a malheureusement des 1* apparition du d§crochage, une zone 
critique oil varie £ortement, et done dgalement, et oil de 

plus, 1' experience donne des valeurs entachees de grosses incertitudes. 
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L* imprecision sur (d$ dans cette zone sera en general indvitable, par- 
ticulierement si I'on s'ocojpe d'un moment. II faudra quand meme essayer 
de le representer au mieux. 

La mesure de doit done etre conduite avec beaucoup de soin, pour 

im nombre elevd d* incidences, garantissant la continuite de 1’ evolution. Une 
solution consiste & utiliser une rampe tres lente d' incidence, qui donne une 
mesure continue. 

La complexity de 1 • aerodynamique ddcrochee, fait que 1 'utilisation d'lme 
telle rampe va donner un rdsultat diffdrent selon le sens dans lequel elle 
est conduite. Sous la forme proposee, le modele ne peut pas rendre compte de 
ce comportement . L’utilisateur devra utiliser une courbe statique moyenne 
qui lui semblera satisfaisante. 

5.2 Moddlisation en configuration non dgcrochge 

II s'agit de rechercher la valeur des coefficients ^ de 

• 

5* A ^ -t- (I A •rV~J & & 

dans tout le domaine non decrochd. 

Cette recherche est faite d partir de la rdponse harmonique de a 

une vibration sinusoidale de petite amplitude. 

Pour les thdories adrodynamiques lindarisdes classiques, cette rdponse 
est inddpendante de 1* incidence moyenne du profil. L'expdrience confirme d 
peu pres ce fait, ce qui nous conduit & utiliser des coefficients r* 

ne ddpendant que du Mach. Cette simplification ne rend d'ailleurs pas toujours 
I'dquation lindaire, puisqu’on est parfois amend a utiliser une function 

(&) non affine f)our le ) . 

Nous avons vu au chapitre prdeddent les rdles des diffdrents paramdtres. 
La quality des rdsultats expdrimentaux enregistrds dans cette zone, la signi- 
fication physique tres prdcise de chaque parametre, font que la recherche 
de ^ peut facilement se conduire manuellement. Nous utilisons 

toutefois I'algorithme ddcrit a 1* annexe A qui automatise cette recherche. 
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5.5 Modilisation en configuration dicrochee 

H C a, ej est une fonction que I’on ne sait pas choisir pour 1' ins- 
tant. Un retard dependant de 1' importance de ^ au passage de 1* incidence 
de decrodiage dewait donner de bons rdsultats. En fait, nous utili- 
sons le retard universel ddfini par Beddoes [8] pour des oscillations 
harmoniques et qui vaut Atm ou simplement A7 = 5 en temps 

reduit (exactement 4,88 avec un dcart moyen de 0,98 sur 142 exemples choisis 
par 1' auteur avec differents profils, frequences, amplitudes de vibrations, 
nombre de Mach) . 


Ce retard est celui de 1* apparition d*un tourbillon au bord d'attaque 
qui fait chuter le moment. La portance ne chute pas tout de suite, mais avec 
un retard moyen de 10,8 en temps reduit, toujours selon Beddoes, d caiise du 
temps de transit de ce tourbillon, mais cet effet est cense etre pris en 
compte par les equations. 

Le choix de la fonction H (&,&) est independant du reste de la mode- 
lisation. 


II nous faut dgalement ddteiminer les valeurs des six coefficients 

A > j , ^ , -t , a » -e des deux Equations. 


Six paramdtres,celd fait beaucoup, mais en se rapportant au chapttre 
precddent, plusieurs simplifications sont possibles. Nous avions conclu que 
la deuxieme equation avait un role preponderant, sauf aux frequences dlevees 
oH Q/© admettait 1' asymptote . Le role de rl en particu- 

lier etait negligeable. 


De plus, il faut considerer la premiere equation comme donnant le Cf 
en incidence non decrochoe, et ce, meme aux incidences supdrieures a 1' angle 
de ddcrochage. II existe des artifices qui retardent le ddcollement : aspira- 
tion de la couche limit e oD fort mouvement de tangage, ou numeriquement , 
resolution des equations adrodynamiques simplifides. On constate alors me 
rdponse qui est dans la continuitd de celle obtenue aux plus basses incidences 
La loi = 2TT^(0en Radians) en particulier semble toujours valable. 


Nous tacherons done de conserver & la premiere equation son comporte- 
ment non ddcrochd, en ne rdservant que quelques corrections sur ces para- 
metres . 
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Nous nous proposerons done s implement de : 

. extrapoler ^ dans la zone ddcrochde. Les points experimentatix ne peuvent 
pas permettre de distinguer un comportement particulier de /I . 

Comme on choisit Co^st ^ 1 * extrapolation est toute faite. 

. corriger Idgerement *r et J dans la zone decollee pour tenir compte de 
tendance d haute frequence que montrerait 1 * experience . 

Si la correction n*appara£t pas sur les points expdrimentaux (comporte- 
ment asymptotique attaint d frequence trop elevee, tout particulierement pour 
^ ), on conserve alors dvidemment les valeurs non ddcrochdes. 

Recherche des coefficients , 3 , e , 


Cherchant a approcher au mieux les points expdrimentaux, on dSfinit 
une fonction erreur dgale a la somme des carrds des erreurs en chaque point 
de mesure, en affectant dventuellement des poids pour chaque frequence, ou 
pour les parties rdelles par rapport aux parties imaginaires. 

Le programme de 1 'Annexe B cherche les valeurs ^ J , c qui minimi- 
sent cette erreur, pour chaque valeur du Mach et de 1' incidence (optimisation 
par la mdthode des gradients ddflechis. 

On pourrait maintenant lisser les , a , e obtenus par des fonction; 
du Mach et de 1' incidence. Mais les lissages que I'on peut faire sont diffi- 
ciles et conduisent a de mauvais resultats du fait de la dispersion des 
valeurs A. , 3. , ^ optimisdes. 

II faut distinguer 1' action des trois variables, en se reportant au 
chapitre precedent . 

Le pic de rdponse de ^ est le caraetdre de la courbe le plus sur, 
mesure a 1' experience. II dvolue peu d'une incidence d I'autre, d'un Mach a 
1' autre. Nous avons dit que sa position dtait ddfinie principalement par 
A. , dgalement par a qui joue lui surtout sur sa hauteur. 

Dans un premier temps, nous nous contenterons done de lisser les 
valeurs de , par une expression -A. . La nouvelle valeur 

peut etre considdrde comme une meilleure approximation de la valeur de 
que celle obtenue par le lissage brut. C'est elle que I'on utilisera ddsor- 


mais. 
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Toujours avec le programme de 1 'Annexe B, on recherche ensuite le couple 
optimum (s,e.) , A. dtant fixe d la valeur que I'on vient de determiner. 

Les differentes valeurs de a sent lissdes par une expression a (<M, . 

Encore une fois, cette nouvelle valeur de j peut etre consider ee comme une 
meilleure approximation de a que celles obtenues par lissage brut, et elle 
sera ddsormais utilisde. 

II ne reste plus que e , dont on recherche dgalement les valeurs 
optiraales lissdes en fonction de et de ^ . 

La modeiisation est ainsi achevde. 

5.4 Remarques importantes concemant le choix des lissages 

II nous faut trouver des expressions en fonction de 1’ incidence, vala- 

bles dans les deux domaines ddcrochd et accrochd. II faut que ces expressions 
ne varient pas dans la zone accrodide, et varient peu aux trds fortes inci- 
dences oil les configurations d’dcoulement se ressemblent. En fait, il appa- 
rait plus simple d'utiliser comme variable au lieu de & . Et comme 

certaines variables ont un comportement trop violent, nous utiliserons 
plutot qui semble la meilleure estimation du ddcollement. 

Ainsi, d correspond I'ordonnde d I'origine des courbes 

</ , yrf , v' , -4 , J , e qui sera la valeur de la variable 

correspondant a I'dcoulement accrochd. 

Les functions -I et ne doivent jamais prendre de valeurs ndgatives. 

Cel& signifierait que le profil h une certaine incidence fixe serait soumis 
a des forces divergentes vers 1* inf ini. Cela ne serait acceptable que pour le 
cas du cylindre oil une force oscillante se ddveloppe spontandment pour cer- 
taines vitesses de vent, mais il faudrait alors introduire des termes non 
lindaires pour empecher cette force d'aller a I'infini par I'dtablissement 
de cycles limites (tout ceci est ddveloppd par Szdchdnyi [3]). 

____ Nous ne pouvons pas isoler le comportement de la premidre dquation dans 
la zone ddcollde, et il a fallu supposer qu'il ne changeait pas de fagon 
importante. De meme, le comportement de la deuxieme equation n’est pas acces- 
sible dans la zone non ddcollde. 

Or, quand I'dcoulement raccroche, 1' excitation de la deuxieme dquation 
s'annule, mais non : la deuxieme dquation continue de rdpondre tant que 

ne s' est pas amorti, ce qui se fait d'une fagon qui ddpend de A, 
de 3 , . 
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II nous £aut done definir Jl et a dans la zone accrochde egaleraent, 
mais les vibrations a petite amplitude ne peuvent alors nous renseigner a leur 
su j et . 

Des experiences spdciales sent necessaires, par example des echelons 
d* incidence de forte amplitude entre la zone ddcrochee et la zone accrochde. j 
La fagon dont le rdgime permanent s’dtablit ddpend de -t et de a . 

Jusqu'5 prdsent, nous avons utilisd les grandes oscillations d'incidence 
a cheval sur 1' incidence de ddcrochage. La fagon dont se fait le recollement 
pendant les boucles ddpend dgalement de -t et de a , mais d’une fagon plus 
complexe. 

Un peu de tatonnement ainsi bien sur que la condition de continuitd de 
■t et de 5 avec les valeurs du ddbut de ddcollement (n^ais ce critere 
n'est pas tres prdcis), mais a conduit a utiliser jusqu’ici 


JL 0.0100 

5 corapris entre 0.15 et 0.20 


5.5 Rdsumd de la mdthode 
. Moddlisation des courbes statiques 

. Moddlisation pour les incidences non ddcrochdes (programme de 
. Moddlisation pour les incidences ddcrochdes 
. recherche des expressions ^ 

. optimisation des triplets d'oii on ddduit 

. optimisation des doublets (^,e) d*oh on ddduit 

. optimisation de e d'oii 1’ on ddduit 


1' Annexe A) 


^CcU,b) 


Les optimisations sont conduites pour chaque Mach et chaque incidence 
par le programme de 1' Annexe B. 

. Vdrification sur des grandes boucles empidtant sur la zone ddcrochde des 
choix des valeurs ^ et A CcM^ . Correction 

dventuelle de ces valeurs, avec la contrainte de rester compatible avec les 
rdsultats expdrimentaux obtenus S petite amplitude. 

Ces grandes boucles permettent aussi de juger de la qualitd gdndrale 
de la moddlisation et de la ndcessitd d’introduire la fonction H L^,B) . 
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6. EXEMPLE DE MODELISATION 

Exemple relatif a la portance du profil OA 209 pour des mouvements de 
tangage autour du quart avant. 

6.1 Definition des variables 



~ variables introduites pour modeliser les couches statiques 
Les autres variables sont definies comme au chapitre 5 page 11 

6.2 Points de mesure 

Ils ont it§ obtenus par le CHAT [4]. 

Le domaine d' application du modele va jusqu’a 20° d' incidence pour des 
Mach compris entre 0. et 0.30, mais son extension jusqu'a Mach 0.4 devrait 
etre possible. 

Des essais a petite anq>litude de vibration ont dte faits a Mach 0.12, 
0.20 , 0.30, pour les incidences moyennes suivantes : 0°, 6°, 9°, 11°, 12°, 
13°, 14°, 15°, 16°, 17° et 18° et mone 20° pour Mach 0.12. 

Les frequences de vibration ont ete : 

pour Mach =0.12 1 .5 - 3 - 5.5 - 10 - 15 - 20 - 30 - 40 Hz 

" Mach = 0.20 3 - 5,5 - 10 - 15 - 20 - 30 - 40 Hz 

" Mach =0.30 4 - 5,5 - 8 - 10 - 15 - 20 - 30 - 40 Hz 

soit des frequences reduites comprises entre 0.045 et 1.23 (M = 0.12) 

0.055 et 0.74 (M = 0.20) 

0.050 et 0.49 (M = 0.30) 
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40 Hz represente une limite physique du montage que l*on ne peut depasser 
sans ditiriorer la quality des mesures. 

Ce jeu de frequences a permis de bien d§£inir toutes les rSponses harmo- 
niques necessaires a la recherche des parametres. 


6.3 Modglisation statique 

Nous avons recherchg les courbes statiques sous la forme : 

pour a decrochage 

*- ^ 7 p°“'' 

ce qui s'ecrit encore : 



Cette famine de courbe s'interprSte ainsi (voir figure 13) : 

. aux basses incidences, suit une relation affine de pente p„ 

. aux tres fortes incidences, il tend vers une asymptote de pente /j , 

(et qui passe par le point d’abscisse dj et d'ordonn6e J 

. enfin, dans la zone decrochee, tend d'une fagon monotone vers son 

asymptote, d'autant plus vite que ^ a une valeur plus negative. 

Les valeurs utilisees qui restituent les courbes expSrimentales sont : 


= 0.03 

p„ = 0.102 / ]/7T^ 

f, = 0 

x: =-0.55 M -«• 0.65 
=-0.3 M -0.43 

=^.21lp, = 12.45 \J'{- Jl''' 

Voir les traces a la figure 14. 

On observe sur les points exp§rimentaux un arrondi au moment ou le 
decrochage se produit tandis que les lissages montrent un decrochage brutal. 
Nous n’ avons pas pris en compte cet arrondi qui amenerait une formulation 
lourde pour un gain ndgligeable. 
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6.4 Modilisation en configuration non decroch§e 

(incidence de decrochage donnee par la formule = 12.45 fuF 

Le programme de 1' Annexe 1 permet de trouver le jeu de parametres 
realisant la meilleure approximation aux moindres carres 

d chaque Mach. 

Un probl^me se pose parfois : les courbes optima ob.tenues passent obli- 
gatoirement par le point d’ordonnde ^ frequence nulle sur la partie 

rdelle. Mais des incertitudes font dans certains cas que ce point ne se situe 
pas bien sur le prolongement des points experimentaux a basse frequence. Les 
lissages sont alors perturbes. Dans ces cas particuliers, nous appliquons le 
programme en imposant le meilleur point de depart possible.^ 

Dans 1* utilisation des parametres ainsi ddterminds, noias utilisons a 
nouveau , car les equations doivent absolument possdder le compor- 

tement statique. Cela decale l§gerement les rdponses harmoniques (voir figure 
15 pour M = 0.12 particulierement) . 

Les parametres retenus qui ont donne les lissages des figures 15,16 et 
17 sont les suivants : 


r\ = 0.20 

= 0.087 

^ = 0.0775 -0.08 M 

Ils seront utilisds sur toute la gamme d’ incidence inferieures a 
1' angle de decrochage. La figure 16 montre I'effet de 1* incidence sur les 
points expdrimentaux. Celui-ci, ndgligd, ne devient important qu*au vois inage 
de 1’ angle de decrochage, dans la zone oil la courbe statique forme I'arrondi 
dont nous n'avons pas tenu compte. 

6.5 ModSlisation en configuration decrbch§e 

Les principes que nous avons deja d6finis sont : 

. fi est prolonge par continuite 

Si les points expdrimentaux montrent des comportements asymp- 
totiques clairs, en deduire les fonctions et 

Sinon, prolonger les valeurs obtenues en non decrochi. 

.recherche de ✓L , puis de a , puis de e 
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Nous utiliserons done McjU>,&) = 0.20. 

Corame les points experimentaux definissent mal la direction asyinpto- 
tique imaginaire, on garde = 0.087. 

II nous £aut maintenant choisir . La, nous avons proedde d'une 

fagon un peu di££erente, de celle etablie qui donne un peu plus de rigueur a 
la recherche du paramStre. Mais il est certain qu*un choix £ait S partir 
d*une estimation de 1* asymptote de la partie rdelle aurait donnS un bon 
resultat . 

7" est trait e comme A y A et e , avec le programme de 1’ Annexe B 

Nous avons soulignd dej& le £ait que I’ordre de la recherche des para' 
metres ne devait pas etre quelconque (voir page 13 ) . Le param^tre a. 
ayant un role plus important que v- , il £audra le rechercher d'abord, pour 
stabiliser les courbes autour de £amilles sures. Nous avons ensuite recherchd 
pour £ixer la hauteur des pics liee & la position de 1’ asymptote 
et dont ddpend beaucoup a et € . Puis nous avons recherchd a CM.&) 

puis e(ctL,e). 

En £ait, nous avons done intercald le choix de r entre A Qt S 
dans I'ordre A , A , e . 

Les rdsultats obtenus sont donnds dans les pages 42 a 48. 

Nous nous attendions a une continuity acceptable en Mach, mais nous 
avons ete contraints de constater une discontinuity entre Mach 0.12 et 
Mach 0.20. Par contre, les rysultats a 0.20 et 0.30 sont su££isamment 
proches pour qu'on puisse utiliser la meme £ormulation pour les deux. 

La di££yrence S Mach =0.12 provient peut-etre d'un e££et du nonibre de 
Reynolds. Toujours est-il qu'un tel e££et est genant pour la dy£inition d'un 
modele. Pour 1' application, nous utiliserons une moddlisation pour les Mach 
in£yrieurs a 0.12, me autre pour ceux supdrieurs a 0.20 et me modyiisation 
interpoiye entre les deux prdeydentes pour les Mach intermydiaires . 

Les rdsultats successi£s obtenus par le programme de lissage sont 
domds aux pages 42, 43 et 44. 

Figure 17 : nous avons lissd en £ait v/r qui est me approximation de la 
position du pic d' amplitude de la rdponse, quand a ne prend pas de trop 
£ortes valeurs. Nous 1' avons recherchy parmi me £amille d 'hyperboles, mais 
des paraboles auraient certainement dty plus simples. 
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L'avantage des hyperboles tenait en leur accroissement plus brut pour 
les tres grandes valeurs de CiC^ qui risquent de se rencontrer en certains 
points de fonctionnement des pales d*h§licoptere. En fait, les autres coeffi- 
cients ont dte lisses par desparaboles et meme une cubique, ce qui annule 
l'avantage escompte. Pour 1’ application, il faudra done faire attention de ne 
pas aller au deld du doraaine lissd. 


avec 



y- oc ^ 

1 

1 pour M ^ 0.12 

1.525- 4.375 M pour 0.12 ^ M ^ 0.20 
0.65 pour 0.20 < 0.40 


Figure 18 : Le lissage de o" n’a pas pose de probl^mes non plus. Une 

famine de droites a suffit. Le point de ddpart de ces droites etait imposd 
a la valeur de non ddcrochde, valeur dSja ddterminee 

7“ = 0.0775 - 0.08 M + T 

avec = -0,19 M^0.12 

/= 1.3875 M - 03565 0.12<M^0.20 

r= -0.079 0.20<M<0.40 

L'examen des tableaux montre bien que I'ordre des choix de , V" , 

<3 et e ne pouvait etre quelconque. Des la premiere dtape, se lisse 

bien mais r prend des valeurs erratiques > alors qu’a la deuxieme 6tape, 
iL se lisse tres correctement. 

Figure 19 : Lissage de a tres ais§ par des paraboles a = t ^ 

Le choix de I'ordonnSe a I’origine vers 0.20 est logique, mais il 
s'avere qu'une tres faible variation de cette ordonnee a I'origine influe 
fortement sur les grandes boucles, au moment oh I'dcoulement raccroche et que 
la reponse de la deuxieme equation s'amortit. Choisir 0.15, 0.20, ou 0.25 
change insensiblement la fonction A. & ) mais peut conduire d des grandes 
boucles doubles les unes des autres. L' experience montre que 1' amort is sement 
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est en general inferieur a 1 ' amort is sement critique qui vaut ici 0.20, d'ou 
le choix de 0.15 pour A 

(I'amortissement critique vaut .3^ = . II est tel que tout retour d 

I'equilibre se fait par oscillations pour a < et sans oscillations 
pour a > Se. ) . 

Les courbes montrent la difference tr§s nette de comportement S 
Mach =0.12. 

On retient 3. = 0.15 - ^ 

avec S = 1.75 pour M^0.12 

f = -16.25 M + 3.70 pour 0.12-e M $ 0.20 

S = 0.45 pour 0.20^ M ^ 0.40 

Les valeurs de e sent tr§s bien aligndes 

£=-2.7 pourM<0.12 

£ = 26.25 M -5.85 pour 0.12^ M ^ 0.20 

£ = - 0.6 pour 0.20^ M ^ 0.40 

II est remarquable que la valeur e = 0 convienne si bien pour 
I'ordonnee 3 I'origine, alors que n' import e quelle valeur aurait etd acceptde. 

Le modele est maintenant compl§tement defini. Les lissages obtenus dan; 
la partie ddcrochde sont donnds aux figures 21 et 22. L* ensemble des rdponses 
hannoniques est montrd par les figures 23 et 24. On y ranarque la discontinuiti 
apportee par le decrochage. Cette discontinuitd aurait 6te remplacde par une 
variation brusque et continue, si on avait tenu compte de I’arrondi des 
courbes statiques. II est remarquable que les points expdrimentaux montrent 
clairement cette variation continue dtalee sur les 2° prdcddant le decrochage 
(voir figure 16, page ). 
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6.6 Les grandes boucles 

Nous avons trace aux pages suivantes, trois series de grandes boucles 
(figures 25.26.27). 


Mach 

0.12 


fryquence 

4 Hz 

; fig. 

Mach 

0.30 


II 

4 Hz 

; fig. 

Mach 

0.30 


II 

8 Hz 

; fig. 


L' amplitude de vibration est de ^ 6° autour d*une incidence moyenne 
comprise entre 0° et 20°. Les boucles autour de 20° d* incidence sortent du 
cadre dans lequel la modelisation a etd faite, mais elles ne conduisent pas a 
un resultat absurde. 

Les calculs ont dte effectuds avec le modele brut, et en tenant 
compte du retard au ddcrochage utilise par Beddoes . La diffirence entre 
les deux jeux de resultat est assez faible, mais elle est certainement dans 
le sens d'une amelioration. Le retard AZ de 5 en temps reduit est assez 
faible. 

Ces resultats devraient etre compares a des essais en soufflerie. Cette 
comparaison aurait dventuellement conduit a une modification des parametres 
a et pour » o . 

Cette comparaison nous a pose quelques problemes qui ne sent pas 
I'objet de ce rapport proprement dit. Elle a done et6 reportde dans 1’ Annexe C 
page . 


7. CONCLUSIONS 

Nous avons expose dans ce rapport une mdthode gdnerale utilisable pour 
la recherche d’un module du comportement a petite amplitude de vibration 
valable sur tout le domaine d*utilisation d*un profil. Les resultats obtenus 
sur un exemple montrent comportement en continuity des coefficients tel 
qu’il valide le choix des Equations utilisdes. 

Les travaux precedents [1 et 2] montrent la possibility pour ce jeu 
d' equations de retrouver le cemportement d grande amplitude, si on se restreint 
aux mouvonents relativement lents, qui forment d'ailleurs les cas d'applicatior 
realistes. 


Nous pensons done que notre modele est rentry dans une phase d’ utili- 


sation. 
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Les travaux devront neamoins se poursuivre par la modelisation 
d*un profil (0A9) sur la gairane complete de Mach et d' incidence rencontree 
sur une pale d*hdlicopt§re, en portance et en moment. Des difficultds peuvent 
apparaitre sur le moment, qui seraient dues h la grande dispersion des points 
experimentaux . 

Des ameliorations sont toujours possibles, qui ne remettraient pas 
en cause les rdsultats jusqu'ici obtenus. 
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ANNEXE A 


Algorithme de recherche des coefficients dans 
le domaine non dScrochd 


Dans le domaine non decroche ^ est solution de I’dquation 

CtJ A C^ = A (J^ -t 7-J 6 & 


A, . ~ 




filers 


^ ^ (Aa-tt) -'AIa 

Crf 9 ^ * 


& 


/ 




On utilise 1’ algorithme de Dat-Meurzec[10] pour optitniser les points expd- 


rimentaux par cette fraction rationnelle 
Notons les points expdrimentaux 


NQIJ 

yot? 


f^l , = 5, = X. .//{ 


Notre famine de parametres reprdsentera d’autant mieux les points 
experimentaux que I'on minimisera 


L Tioi; . 


erreur= ^ ^ j 


ou 7^ est un poids positif 


et 


designe la norme ^ ^ ^ 
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Application de cet algorithme a la deuxieme Equation 

Le meme raisonnement nous conduit au systeme suivant : 

Zb. n li\ = Z B 

* ^ [el ^ 

1 pour la premiere iteration 

(V - 1 ) M. o. <a.et-^ etant les valours 

trouvdes a 1’ iteration precedente. 


avec ^ = 



Cet algorithme ne permet pas de prendre V en compte comme nous 
I'avons fait dans 1’ example traitd, puisque la rdponse n’est plus sous 
la forme d'lme fraction rationnelle mais d'une somme de deux fractions ration- 
nelles. Si on ramene cette somme & une seule fraction, les coefficients 
inconnus interviendront pas leurs produits d'une fagon que I'on ne peut 
traiter. 

L' application brute de cet algorithme a toujours conduit a des 
calculs divergents CcC , 
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ANNEXE B 


Optimisation des parametres par la minimisation 
d’une fonction de plusieurs variables 


Notre probleme consiste & lisser une serie de points experimentainc 
par une fonction de transfert qui depend 

de r> parametres (en 1’ occurrence ^ , T , 3. , e ) 

Pour cela, nous d6£inissons une erreur 


etant un poids eventuel qui depend de la frequence rdduite, 
et un autre poids destine S donner plus d' importance a la partie 

r§elle ou a la partie imaginaire. 


La mdthode consistant & rechercher les racines du systeme 


? { Ifc) 

Jfi 


= 0 


-V« 


conduit a des syst§mes impossibles a resoudre directement. 

Nous avons done recherche le minimum de 1’ erreur par une mithode 
de gradients construite autour des formules 22 pages 2.14 de la [6] 


1 . Algorithme 

1 . Choix d*un point de depart Tc pour 


t = 1. 


Nous avons suffisamment 1’ habitude des systemes etudi§s pour con- 
naitre des points de depart valables. 


Pour une portance, nous utilisons : 


= 0.04 

^ = 0.1 

d, = 0.2 

e, = 0 
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2. Calcul du gradient Cc de frr , analytiquement . 

3. Calcul de la direction de recherche Tc qui n'est pas celle de plus grande 

pente , mais une direction amelioree qui tient compte des proprid- 

tds de courbure de la surface ^ = Brr C^i, / '^k) , courbure 

deduite des pas precedents : 


~Tcs -Gc + 


; 

si t = 1 , on fait 71- = - 


HSc.,t 




Enfin, pour une meilleure consequence, nous reinitialisons 'Tc ^ 
toutes les 7 iterations. 

4. Recherche de tel que I'erreur au point TV- + ^ *1’^ , 

notde FU) soit minimum. 

On recherche done ^^1.-0 soit 


4^ (r^,- 

J / 




= 0 


( (jc)j' ddsigne la composante de 7T-J 

V 

4 a. Premidre approximation de • 


P^tet 

utiUse) 


On choisit une valeur ^ telle que soit la plus grande puissance 

de 2 infdrieure S , valeur que l*on di vis era ensuite par 2 pour 

avoir une meilleure sdeuritd du calcul (pour dviter de se trouver trop loin 
au deia du minimum recherchd, dans vine zone oh frr risque de retrouver un 
comporteraent complexe) : 


/WT 


4 






INT = function partie entiere. 
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4 b. Recherche d'un encadrement de la solution : 

A partir de la propriete ^ 

dL 

- si on divise ^ par 2 

jusqu'a ce que ‘“'/4c <0 


- si dL f£^£.') on double ^ 

dt ' 

jusqu’a ce que df/4l -^0 

Dans les deux cas, on retient les deux demieres valeurs de 
rencontrdes comme homes de I'intervalle ["/j, cherche. 

Le deuxione type de boucle (on double ) risque de ne 

jamais s'arreter si le minimum se trouve sur une direction asymptotique^ 

(i-Dans le cas de notre modelisation, certaines optimisations correspon- 
dent effectivement a une valeur infinie. II n’y a alors rien a faire. 


4 c. Reduction de I'intervalle en 

plus petit, par dichotomic, en se basant sur le signe de 
choisir I'intervalle inferieur encadrant . 

4 d. Interpolation lingaire sur pour trouver , 

approximation de 



5. Calcul du nouveau point £p.~n. 

et retour en 2, sau£ si la distance est suffisamment 

faible, auquel cas on considere la convergence terminge ( T;T*c*i 10 ) 


if, pour 

d£ 


notre meilleure 
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2. Rgsultats 

Notre test d' arret de convergence est tres s6v§re ( ^ 10 

mais il n'y a pas de risque de divergence au voisinage du minimum, done autant 
employer un test difficile qui assure mieux la proximity de la solution 
cherch^e. D' autant que nos calculs sont effectues sur une machine de bureau 
et sont gratuits 


Nonibre d’ iterations effectuges en moyenne : 


4 degr§s de libertd 


3 

2 

1 


tt 

11 

M 


erreur faible 


i 20,25 

i 10,12 
- 8,10 
2 


erreur forte 
jusqu'§ 80 
jusqu’a 25 
jusqu'a 15 
2 


Ces rgsultats ne sont que des indications relatives auccas traites. 


3. Am§liorations possibles 

Les formules utilisees en 3 ameliorent nettement la consequence 
du systeme. D'autres formules proposant des directions de recherche meilleures, 
peuvent etre utilis§es (voir en particulier [7]. 

A no ter qu’ a la[6]» page 2.15, figurent aussi les formules 23 qui 
ont etg utilisees ggalement mais qui ont donn6 un. mauvais resultat (il y 
figure probablement une erreur) 





RT 12/1841 RY 080R 


Fage 65 


ANNEXE C 


Comparaison du modele avec 1' experience pour I'exemple qui 


a servi d' illustration 


Le trace des grandes boucles calculdes (figures 25, 26, 27) montre 
un ddsaccord avec les grandes boucles mesurdes (figures 30, 31 et 32, lignes 
pointilldes) . Un rdsultat analogue avait dtd obtenu par Mile Falchero [9] . 


1 . Probleme de la qualite des essais en soufflerie 

La cause principale de ce ddsaccord semble provenir du fait que les 
essais grande amplitude et les essais petite amplitude, proviennent de detix 
campagnes distinctes, et ont dtd rdalisdes sur deux maquettes diffdrentes du 
profil 0A9. 

La premibre campagne [5] port ait sur la mesure de la rdponse du 
profil a divers mouvements harraoniques, grande amplitude, mais aussi petite 
amplitude (a 20 Hz , 30 Hz , 40 Hz) . La maquette avait ensuite eti ddtruite pour 
rdcupdrer les capteurs pression. Le processus de modeiisation demandant 
beaucoiQ) plus de point de frequence, une deuxidme campagne [4] a eu lieu, qui 
a conduit au modele du chapitre 6. 

Une certaine dispersion apparait bien sQr entre les deux campagnes, 
chaque fois que les memes quantites ont ete mesurdes : 

. figure 28 : les courbes statiques sent un peu diffdrentes, avec un ddcro- 
chage qui apparait un peu plus tot a la deuxieme campagne et un ddcroche 
moins fort, line troisieme courbe statique, obtenue S Modane est dgalanent 
reprdsentee, qui donne un rdsultat encore different, avec un angle de ddcro- 
chage particulieranent dlevd. 

. les mesures S petite amplitude de vibration communes (frequences de 20,30 
et 40 Hz) montrent des variations dgalement. On peut voir en particulier, 
grace aux tableaux de la figure 29 que le maximum de la partie rdelle qui se 
situe entre 15 et 20 Hz lors de la deuxieme campagne serait plutot vers 30 Hz 
a la premiere campagne, ce qui correspondrait a des valeurs de un peu plu 
fortes . 
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D'me maniere g^nSrale, les essais en soufflerie tendent a se dispersei 
au mcsment du decrochage et dans la zone ddcrochde. II su££it de peu de chose 
par exemple pour modi£ier un angle de dicrochage, ce peu de chose pouvant 
etre une irrdgularitd de pro£il, ou un e££et tridimens ionnel. II est possible 
que les petites vibrations soient des mesures plus sures que les courbes 
statiques, mais ce n'est pas evident. 

Tout ceci est assez ennuyeux car les di££drences que I'on note sent 
signi£icatives dans le module. L’utilisation des courbes statiques de la 
premiere campagne, coup lee avec les coe££icients ddduits de.la dexjxieme cam- 
pagne amdliore grandement nos boucles (voir £igures 30,31 et 32). 

Le modele tel qu'on l*a dd£ini schematise correctement le comportement 
d'une certaine maquette dans sa sou££lerie. Les de£auts du pro£il de la 
sou££lerie, des traitements des mesures (1 Vage du sou££leur a-t-on 1 'habitude 
de dire . . ) sont aveuglement pris en compte dans nos equations , qui rendent 
compte du comportement grande amplitude du pro£il dans les memes conditions 
de mesure. Ces conditions sont choisies pour etre les plus proches possibles 
de la rdalite. Mais les sou££leries sont des moyens d' essais complexes d 
maitriser. Bien evidemment, ces probl^mes sont ceux de tout module s’appuyant 
sur des essais sou££lerie. 

2. Application du modele aux grandes boucles de la premiere campagne CHAT 

Pour cette application, il est necessaire de tenir compte de certains 
resultats obtenias a la premiere campagne, au moins des courbes statiques que 
le modele doit en particulier prevoir. 

Nous avons done recherche pour ces courbes statiques des £ormules 
analogues a celles ddja utilisdes : 




^= 0.05 
X ^ 0.102 
0 




4/3 M + 0.7 
^^- 0.5 

^^16° - 10 M 


ec&j 
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C’est cette unique correction qui a donne les figures 30, 31 et 32. 
Les resultats sont maintenant acceptables et auraient certainement meme ete 
un peu meilleurs si tous les essais avaient §te faits a la suite les ims 
des autres, tout particulierement aux tres fortes incidences. 

Nous avons regardd plus complgtement I’effet des corrections sur une 
boucle a cheval sur 1* incidence de ddcrochage. Se reporter § la figure 33. 


Les conditions sont 


Incidence 

: 11° + 

Mach 

: 0.30 

Frequence 

: 4 Hz 


lere boucle : application brute du modele obtenu d partir de la deuxieme 

campagne. Le resultat n’est pas bon. 

2eme boucle : prise en compte du retard de Beddoes qui prolonge la partie 

droite de la boucle, mais n’amdliore pas globalement . 

3eme boucle : utilisation du ^ statique de la premiere campagne. La 

boucle est maintenant bonne. 

4eme boucle : modification de pour tenir compte de la remarque apportde 

par les tableaux de la figure 29 : 


0.1 + 0.65 bCy - 1 + 


0.65 + 1 


devient 


0.15 + 0.65 - 1 + 


0.65 - 1 


Pour ne pas changer le comportement de 1 ' aSrodynamique , on augmente 
d dgalement, de mani^re a garder la quantity a/-t sensiblement constante 
(ce qui maintient un am ort issement un peu infdrieur a I'amortissement critique 
en ^/^ }) 

<3 =0.15 - 0.45 devient a = 0.25 - 0.45 
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La reponse de la deuxieme dquation est m peu amelioree dans la zone 
od 1 * dcoulement vient de se raccrocher 

Seme boucle : pour avoir ime boucle parfaite, on augmente un peu le retard 

A? en prenant 8 au lieu de la valeur 5 proposee par 
Beddoes . 

6eme boucle : comme les rdsultats expdriraentaux dont nous disposons sont 

coupds & partir de la 4eme harmonique comprise, on fait subir 
le mime sort a notre riponse. La modification apportie par 
ce filtrage est assez importante. 

Le resultat final est excellent, sans que I'on puisse conclure bien 
sur que ces manipulations soient nicessaires. Elies ditirioreraient peut-etre 
d'autres boucles. 

Sur tous ces risultats, on remarque a chaque fois que nos boucles sont 
toujours dicalees vers le bas. Est-ce un probleme de statique encore ? 

ou un effet de la seconde equation n’apparaissant pas en realite ? Aux 
frequences reduites exagirees (n’oublions pas les limitations en frequence 
amplitude du modele), les boucles calculies sont tres fortement decentries vers 
le bas. Cela se produit d’ailleurs parfois pour M = 0.12 a une friquence de 
8 Hz ou la frequence reduite de 0.25, ilevee du fait du I^ch tres faible, 
est trop grande pour une amplitude de 6°. 
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